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Abstract 

We present a determination of the neutrino mass matrix which holds for values 
of the neutrinoless double beta decay effective mass m ee larger than the neutrino 
mass differences. We find eight possible solutions and discuss for each one the cor- 
responding neutrino mass eigenvalues and zero texture. A minimal structure of the 
perturbations to add to these zero textures to recover the full mass matrix is also 
determined. Implications for neutrino hot dark matter are discussed for each solution. 



1) Introduction : The neutrino mass matrix in the three neutrino framework has been 
severely constrained over the last few years by atmospheric |l] and solar || neutrino exper- 
iments, and by recent positive indications from laboratory experiments |3|, In the basis 
in which the charged lepton mass matrix is diagonal, the neutrino mass matrix contains 
all the information about the mixing schemes in the leptonic sector. The mixing angle 
required for the atmospheric neutrinos appears to be maximal or very close to maximal. 
This is one of the main constraints on the neutrino mass matrix, when combined with the 
mass squared difference required for atmospheric neutrinos 5m? atm , which is in the range 
(1.5 — 4.8) ■ 10 _3 ||. There are a few solutions to the solar neutrinos, which give mass 
squared differences 5m 2 sol and mixing angle 9 with u e which can lie in various ranges (see 
e.g. Ref. f|). The reactor experiment CHOOZ gives in addition a strong upper limit on 
the U e 3 element of the neutrino mixing matrix 0. Besides these five constraints on the 
neutrino mass matrix, a sixth one would be given by the measurement of the effective 
mass m ee in neutrinoless double beta decay. This parameter would provide a scale to the 
neutrino mass matrix. Very recently a positive indication of neutrinoless double has been 
reported with m ee in the following range ||: 

m ee = (0.05 - 0.86) eV (at 95% c.l.) (1) 

and with a best fit for m ee = 0.39 eV. In the following we study what could be the 
consequences for the neutrino mass matrix of a value of m ee within the range of Eq. (fj). 
By considering this range we study more generally what would be the effect on the mass 
matrix of an effective mass m ee determined with a value greater than the neutrino mass 
differences and smaller than the dark matter constraints (see below). This range covers 
also the values of m ee which could be further tested experimentally in a relatively near 
future |J. Since the scale m ee is then higher than the mass squared differences required by 
the neutrino oscillation results, this implies a degenerate or partially degenerate solution 
pr^| . The hierarchical solution is not allowed by a value of m ee within the range of Eq. (fj) . 
An inverse hierarchy solution is still possible but only for m ee close to the lower edge of 
this range. 

Considering the fact that the atmospheric mixing is close to maximal, several texture 
mass matrices have been constructed Jll[|, including in particular the ones where the mass 



squared difference Sm 2 tm is considered to be the dominant entry in the mass matrix and 
where the solar neutrino solutions could be added as perturbations. These textures are 
very useful to construct models, where the zeroes could have their origin in some symmetry 
of the model, so that the perturbations come from the symmetry breaking. For a value of 
m ee within the range of Eq. (ffl), 5m 2 atm would likely becomes a perturbation in the neutrino 
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mass matrix with respect to the scale m ee . This would be an important information for 
the determination of the texture of the mass matrix. Moreover, given a value of m ee in 
the range of Eq. (|I]) the neutrino contribution to the hot dark matter would be significant. 
With these assumptions we now attempt to look into this problem. Instead of starting the 
analysis from the possible textures, we start with the experimental inputs and then try to 
determine the neutrino mass matrices and in the process get the possible textures. 

In the following we shall work in the basis in which the charged lepton mass matrix is 
diagonal. The physical states \v a >, (a = e,/z,r), which enter in the charged current, are 
related to the mass eigenstates >, (with masses rrii, z = 1,2,3), by the mixing matrix 

\u a >= U ai \Vi > (2) 

and the effective mass entering in the neutrinoless double beta decay is 

m ee = (M u ) ee = (M„) u = U ei m i- ( 3 ) 

i 

{My)u is the (11) element of the neutrino mass matrix in the flavour basis. 

In this article we do not consider the imaginary part of the mass matrix assuming that 
there is no CP violation, or very small CP violation. In this limit, we can parametrize 
the mass matrix in the flavor basis as a function of the three real mass eigenvalues and 
the three angles coming from the mixing matrix U ai . For simplicity and also as a first 
step of a more involved analysis we consider in addition only the case with exact maximal 
mixing between the and v T (motivated by the best fit value from SuperKamiokande) 
and we consider ?7 e 3 = (motivated by the result of the CHOOZ experiment). In this 
framework we determine what are the effects on these six mass matrix parameters of the 
five experimental constraints mentioned above taking a given value of m ee . We find that 
there are eight different solutions for the mass matrix, each one being determined by these 
six constraints. For each solution we display the corresponding zero texture and show how 
the perturbations arise around these zero textures to give the full mass matrix. We also 
discuss the consequences for hot dark matter for the different solutions. 

2) The neutrino mass matrix solutions : In a general way the diagonal neutrino mass 
matrix can be written as: 
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here 02,3 are the Majorana phases, which introduces CP violation in lepton number vio- 
lating processes. The mixing matrix U a i may be parametrized as 



U 




ClC 3 

- C!S 2 s 3 e 



-S1C3 
l<t> cic 2 + sis 2 s 3 e l<i> 




(5) 



cic 2 s 3 e^ cis 2 - s 2 c 2 s 3 e 

where, = sin#i and q = cos^. This depends on 3 mixing angles 9i and one CP- violating 
phase 0. Thus there are 9 parameters in the general mass matrix in the flavor basis |T2 . 



The atmospheric neutrinos require an almost maximal mixing between the flavour states 
and v T . A very little mixing with the v e is not ruled out, but it is constrained by the 
CHOOZ result. Considering this constraint and the fact that the central value of the mixing 
angle for atmospheric neutrinos is maximal, as required by the SuperKamiokande result 
and supported by the K2K preliminary data, we assume only maximal mixing s 2 = 1/ v2 
between the flavour states and v T . Although the CHOOZ constraint allows a small U e3 , 
for construction of a simple form of a mass matrix we assume it to be vanishing, S3 = 0. 
In this limit of U e3 = the CP-violating effect from the phase in U vanishes and we can 
neglect its effect. We can then write down the mixing matrix in terms of the mixing angle 
9 = 81 required for the solar neutrinos solutions as 



U 



( c -s \ 
\ s/V2 c/V2 l/y/2 J 



(6) 



with c = cos 9, s = sin 9. 

This leads to the following mass matrix in the u e , u^, v T flavor basis: 



M,, 



I c 2 m 1 + s 2 m 2 ^(mi - m 2 ) ^(mi— m 2 ) \ 

m 2 ) - 2 ^ 

m \ s 2 rni+c 2 m,2—m3 s 2 mi+c 2 m2+m:j 

2 2 / 



V T2^ 



(7) 



In the following we disregard the effect of the Majorana phases 02,3 of the diagonal mass 
matrix. These CP violation does not enter in any neutrino oscillation experiments. We will 
consider only real (positive or negative) mass eigenvalues. By considering both positive 
and negative mass eigenvalues we consider the extreme cases, and the case with non-zero 
phases are intermediate cases between these two extremes. The results with phases can 
be obtained as a simple generalization of what we present here. Under these assumptions, 
to cover all physical distinct possibilities without loss of generality, it is sufficient to take 
m ee positive and \m 2 \ > |mi| by definition and to consider the range 0° < 9 < 90°. Then 
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from the value of 5m 2 tm , Sm 2 ol , 8, and m ee we can determine the mass matrix. There are 
8 solutions. To find them we have first to solve the equations: 

m\-m\ = 8m 2 sol (8) 

m ee = mic 2 + m 2 s 2 (9) 

This gives two sets of solutions: 

m? = m ee -^^ (10) 



2 m 



c 5m 



ee 



2 Krr-t^t 2 

sol 



and 



TO? = m ee + ^ 11 

2 2 to rr 1 ' 



m? = ^ + ^^ (12) 



cos 28 2 m, 



ee 



mi = -J^-l S <* (13) 
cos 28 2 m ee 

where one can safely keep only the first order in 5m 2 sol /m 2 e as we did here. Then we can 
determine the mass 777,3 from the equation: 



m 



777? + 777 2 



2 

From this we find four solutions with the set mf 2 : 



■,-^^±6ml tm (14) 



mf uM = Jm 2 ee ±5m 2 atm (15) 



771 



3 3 ' A4 = -\m 2 ee ±6m 2 atm (16) 



and four solutions with the set 777f 2 : 
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Here in 7773 we kept only the correction in 5m 2 atm /m 2 ee and neglected the tiny correction 
from Sm 2 ol . 

3) Neutrino mass numerical values : Putting Eqs. fllOD- (|l8D in Eq. (|7|), we can directly 
build the eight possible matrices and substitute the values of 777 ee , 5m 2 atm , 8m 2 ol and 8 to 
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mi = 0.390 




mi = 0.848 




m 2 = 0.390 




m 2 = -0.848 




m 3 = 0.394 


B x 


m 3 = -0.850 




0.050 < mi < 0.86 




0.0750 < m x < 7.74 




0.050 < m 2 < 0.86 




-7.74 < m 2 < -0.0750 




0.0632 < m 3 < 0.8628 




-7.74 < m 3 < -0.0845 




mi = 0.390 




mi = 0.848 




m 2 = 0.390 




m 2 = -0.848 


A 2 


m 3 = 0.386 




m 3 = -0.846 




0.049 < mi < 0.860 




0.0750 < mi < 7.74 




0.050 < m 2 < 0.860 




-7.74 < m 2 < -0.075 




< m 3 < 0.859 




-7.74 < m 3 < -0.029 



Table 1: Central values and allowed ranges of the mass eigenvalues mi j2) 3 (in eV). For the 
solutions A 3 , A4, B 3 , £? 4 , the mass eigenvalues are the same as that of the solutions Ax, A 2 , 
A 3 , A^, respectively, except for the sign of m 3 which has to be reversed. We put several 
digits just to show the effect of Sm 2 tm which is different for the various solutions. 

obtain the possible matrices. To illustrate the differences between the various solutions, 
we give in Table 1 the values of the masses mi j2j 3 obtained in the eight cases, taking, 
as a representative value of m ee , the central value from ||, m ee = 0.39 eV, and taking 
for the other parameters the following central values in the case of the LMA solution: 
Sm 2 atm = 2.9 • 10~ 3 eV 2 [|, 6m 2 sol = 3.7 • 10" 5 eV 2 g, tan 2 # = 0.37 @. The uncertainties 
on these masses due to the experimental uncertainties on these four parameters are also 
given in Table 1. To get them we took the following ranges: 5m 2 atm = (1.5 — 4.8) ■ 10~ 3 eV 2 
H|, 5m 2 sol = (0.2 - 3) • 10- 4 eV 2 §, tan 2 6 = 0.2 - 0.8 [|, m ee = (0.05 - 0.86) eV g. The 
uncertainties on U e3 and the atmospheric angle also induce an error, which, however, is 
moderate and which we didn't consider. We will discuss these central values and ranges in 
sections 5 and 6. 

For the SMA solution, although this solution is now almost ruled out 0, the same 
calculation can be done and for completeness we give in Table 2 the values of 771^2,3 we 
obtain considering the following central values from Ref. 0: 5m 2 ol = 5.2 ■ 10~ 6 eV 2 and 
tan 2 6 = 1.8 • 10~ 3 . We don't give here any range because the central values for the SMA 



6 



A 1 


mi = 0.390 
m 2 = 0.390 
m 3 = 0.394 




mi = 0.391 
m 2 = -0.391 
m 3 = -0.395 


A 2 


mi = 0.390 
m 2 = 0.390 
m 3 = 0.386 


B 2 


mi = 0.391 
m 2 = -0.391 
m 3 = -0.388 



Table 2: Central values of the mass eigenvalues mj 2 3 (in eV) for the SMA solution. For 
the solutions A 3 , A4, B 3 , B4, the mass eigenvalues are the same than for the solutions Ai, 
A 2 , A 3 , A^ respectively except for the sign of m 3 which has to be reversed. 



solution have been obtained in Ref. |3| beyond the 3-cr level. The results for the LOW 
and VAC solar neutrino solutions can be obtained in the same way but we don't present it 
because they are in general difficult to accommodate with a value of m ee larger than the 
neutrino mass differences [pxjl . See however counter-examples in Ref. [fL3| . 

4) Zero textures and perturbations: To understand better the structure of these ma- 
trices and also for the purpose of model building it is interesting to write these matrices 
in the form of textures zero plus the perturbations. This will also allow us to compare our 
results with the textures predicted from different considerations ||11|| . In the following we 
shall consider the degenerate solution assuming that m ee is larger than the neutrino mass 
differences. For the inverse hierarchy solution we will just give a few comments in sections 
5 and 6. 

In the case where the mass eigenvalues are almost degenerate with the mass squared 
differences to be small compared to the eigenvalues, i.e. |7Tii j2 ,3 1 3> ^ m atm ^ ^ m 1oii the 
parameters 5m 2 atm and 5m 2 sol can be considered as perturbations with respect to m ee ~ 
1^1,2,3 1 • I n this case we get the following zero textures and perturbations for the different 
solutions: 
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(19) 
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M A , 



m e 



+ m e 



( 1 \ 

1 ±m e 

V° 1 / 
( 

—\/2sin20 
\ — v/2sin20 



^sol 



/o 







-1 
+1 





+1 
-1 
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— \/2sin20 — \/2sin20 \ 
2 cos 20 2 cos 20 
2 cos 20 2 cos 20 y 



(20) 



M, 



Bi,B2 



cos 20 



( cos 20 

sin 20 
sin 26* 



sin 26 
V2 

1 _ cos 20 

2 2 

1 cos 26* 



V2 



sin 20 
V2 

1 _ cos 20 

2 2 
_1 _ cos 20 

2 2 



m ee COS 2$ 



^atm 



+ m ( 



2 



-1 +1 
+1 -1 

v^sin 29 v^sin 29 

V^sin 29 -2 cos 20 -2 cos 29 

\ >/2sin20 -2 cos 20 -2 cos 20 
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(21) 



M, 



B$,B4 



cos 20 



/ cos 29 

sin 26* 
sin 26* 



V 



sin 20 

V2 
_ cos 20 
2 

cos 26* 



sin 29 
V2 

1 _ cos 20 
'2 2 
L cos 20 



/ 



m ee cos 20 



^atm 



+ m. 



£-sol 



( 



( \ 

1 -1 

V o -i i / 

>/2sin20 >/2sin20 \ 



\/2sin20 -2 cos 20 -2 cos 20 
^ sin 20 -2 cos 20 -2 cos 29 J 



(22) 



Here £ SD 2 = 5m 2 ol /m 2 e and e atm = 5m 2 atm /m 2 e . The corresponding mass eigenvalues are 
given in Table 3. Note that the various solutions differ from each other by their zero texture 
and by prefactors in the perturbations. Except for these different prefactors the various 
perturbations have the same matrix structure for all the mass matrix solutions. 



5) The ^1,2,3,4 solutions: First let us discuss the solutions A 1 to A A : 
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M (0) 


"7.1,2,3 


A lt2 


m ee | 


^10 0^ 

1 

V° 1 J 
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mi = m ee (l - s 2 ^-) 
m 2 = TO ee (l + c 2 ^-) 

TO 3 = TO ee (l ± 






( 1 ^ 
1 

Vo 1 y 
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TOx = TO ee (l - S 2 ^-) 
TO 2 = TO ee (l + C 2 ^-) 

m 3 = -TO ee (l ± 




to j 


I cos 26* 

sin 20 1 
V2 2 
sin 20 1 
V V2 2 


sin 20 sin 20 \ 
V2 v/2 \ 
cos 20 1 cos 20 
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cos 20 1 cos 20 1 
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mi = to(1 + s 2 008 20^-) 
m 2 = — to(1 + c 2 cos 28^f-) 
m 3 — — to(1 ± cos 2 29^r-) 


-83,4 


fh 
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cos 20 1 
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cos 20 1 
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toi = to(1 + s 2 cos 20 ^-) 
to 2 = — to(1 + c 2 cos 29^^) 
TO3 = to(1±cos 2 26>^f^) 



Table 3: Zero textures and neutrino mass eigenvalues at first order in e atrn and e so i corre- 
sponding to the solutions ^1,2,3,4 and to -81,2,3,4- 



• Mass eigenvalues: For these solutions the absolute values of the masses of the 3 
eigenstates coincide with m ee modulo small variations (proportional to Sm 2 tm and 
5m? sol ). This explains why in Table 1 the upper limit of the neutrino mass ranges 
reproduce to a good approximation the upper limit on m ee . For low values of m ee 
(around 0.05 eV in Eq. ([I])), the neutrino masses 777.1,2 are still given to a good 
approximation by m ee but the |m 3 | value can differ substantially from the one of 
m ee . This is due to the fact that in this case 5m 2 tm is not small anymore with 
respect to m ee and |m 3 | could be as low as eV for the solutions A 2 ,4 (where 5m 2 atm 
and m ee come with opposite signs in 777.3). This later case corresponds to the inverse 
hierarchy scenario. 

• Dark matter . For a value of m ee such as the central value m ee = 0.39 eV || or 
above, the solutions ^1,2,3,4 could be relevant for dark matter. \m\\ = \m 2 \ = {m^l = 
m ee = 0.39 eV (or the upper value 0.86 eV in Eq. ([[])) corresponds, in units of the 
critical density, to a value of fl u h 2 ~ Yl |«7j|/(91.5eV) around ~ 1.5% (~ 3%) with 
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h the dimensionless Hubble constant. Note that the value Q u h 2 ~ 3% turns out to 
be about the upper bound allowed by fits from observed galaxy clustering and the 
Lyman a forest (see e.g. |L5|]-[|I8] and the discussion in [|Kj). 

Textures and perturbations in the degenerate case: For the cases Ai t z and A 3A we 
recognize the textures CO and C3 of Altarelli-Feruglio, respectively HTi| (naturally 



leading to both simple maximal mixing or double maximal mixing). From Eqs. (|7|)- 
( |i~8|) one can observe also how the perturbations arise with respect to these zero 
textures to give the full mass matrix. Note that the structure of the perturbations 
due to 8m 2 atm is quite simple, which is interesting for model building. Note also that 
in case v4 lj2 it is this perturbation which is responsible for the maximal atmospheric 
mixing. In case A3 4 the maximal mixing is already present for Sm 2 tm = 0. The 
structure of the perturbation from Sm 2 ol is more complicated but remains rather 
simple especially if 9 takes a specific value (i.e. for example for 9 = or 9 = it /A). In 
all cases A\ to A4 the solar mixing is due to the perturbation proportional to Sm 2 ol . 
The solar angle even if maximal is a perturbation because all the terms depending 
on it are multiplied by Sm 2 ol . There is no solar mixing if 8m 2 ol vanishes. 

The inverse hierarchy case: If m ee turns out to be not much larger than the lower 



edge of the experimental range in Eq. ([!]), i.e. m ee = 0.05 eV, we could have an 
inverse hierarchy where |m 3 | << |m 12 | which means, since m 2 ee » Sm 2 ol , that 
1 777.3 1 << I w ee |- Such a situation requires a cancellation in 777.3, Eqs. (|T5|) and flT5|), 
of the m ee term with the term proportional to 5m 2 tm . This is possible only with 
solutions A 2> 4 where Sm 2 tm comes with an opposite sign with respect to the term in 
m 2 e in 777,3. For these solutions this requires m 2 ee ~ 5m 2 atm . For the central values 
5m 2 atm = 2.9 • 10~ 3 eV 2 and tan 2 9 = 0.37 @ this requires m ee ~ 0.054 eV. Varying 
5m 2 tm within the range (1.5 — 4.8) • 10~ 3 eV 2 this requires: 

4 • 10~ 2 eV < m ee < 7 ■ 10~ 2 eV (for solutions A 2A ) (23) 

This range is independent of the solar neutrino solution. Note also that we could 
develop in this case the mass matrix considering and a/ 8m 2 sol as perturbations 
with respect to m ee and a/ 5m 2 atrn . In this case the zero texture we obtain is 

/ 1 \ 

1/2 1/2 , (24) 

\0 1/2 1/2 J 

that is to say the texture B2 of Ref. |l ij which can lead to both single and double 
maximal mixing. 



10 



6) The -81,2,3,4 solutions: Solutions B\ to B± are very different from A\ to A<± and more 
complicated because mi and m 2 here have opposite sign (c.f. Table 3): 



Neutrino mass eigenvalues: Here in the limit 5m 2 atm = 5m 2 ol = the absolute values 
of the neutrino masses are not m ee (as with solutions A\ to A4) but turn out to 
be in fact fh = m ee / cos 29. This explains the large value of the neutrino masses 
obtained in Table 1 with these solutions. Taking the central value tan 2 # ~ 0.39 || 
for the LMA solution we then get for the neutrino masses m ~ 2.2 m ee . For a value 
of m ee = 0.39 eV this gives fh ~ 0.85 eV. Taking the upper limit tan 2 9 ~ 0.8 |J we 
get fh ~ 9m ee which, for m ee = 0.86 eV in Eq. ([IJ), could result in a neutrino mass 
as large as 8 eV. However, a neutrino mass greater than 2.8 eV is excluded by the 



tritium beta decay experiment ||19|| . 



Dark matter . These values fh = 0.85 eV and fh = 8 eV correspond to Q u h 2 ~ 3%. 
and to Q u h 2 ~ 25% respectively. From Table 1 we then get the following range for 
solutions -81,2,3,4: 

0.001 < tt v h 2 < 25% (25) 

Note that, as already mentioned above, a value of Q u h 2 larger than ~ 3% which 
corresponds to ~ 3 eV is disfavoured by the results of galaxy clustering and 



Lyman a forest experiments |15[-JT8[. In case of the solutions -81,2,3,4, to have Vt v h 



2 



smaller than ~ 3% we need that the following constraint is fulfilled: 

m = ^-<leV (26) 
cos 26 v ; 

For example this condition for m ee = 0.86 eV in Eq. (P requires tan 2 9 < 0.08, which 
can be compared with the lower limit of the LMA solution which is tan 2 9 ~ 0.2 
[^]. This means that a value of m ee as large as 0.86 eV with the solutions -81,2,3,4 is 
disfavoured by dark matter fits. To obtain fh < 1 eV with for example tan 2 29 = 
0.2, 0.39, 0.8 we need m ee < 0.7 eV, 0.4 eV, 0.1 eV respectively. Would we allow 
fl u h 2 to be as large as 10%, these bounds would be multiplied obviously by a factor 
3. Note that for a ACHDM model with both hot and cold dark matter as well as a 
cosmological constant a value of fh < 0.5 eV is favoured [^tj. This already favours a 
lower value for the neutrinoless double beta decay in this scheme. Furthermore, the 
CMB studies with MAP and PLANCK and the Sloan Digital Sky Survey experiment 
(see e.g. |Tj]]-[E2]]]) will be able to probe the mass range fh ~ 0.25 — 0.5 eV. If the 



contribution to the hot dark matter from neutrinos is constrained to be less than 
fh < 0.25 eV, then this will mean that this class of solutions will be further restricted 
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and a lower value of the neutrinoless double beta decay will be favoured. Note also 
that the fact that the LMA solution in the degenerate case could be relevant for 
dark matter is well known Jbl. Here we emphasize that, taking all the experimental 



constraints on U e3 , 5m 2 atm1 5m 2 olar and taking m ee within the range of Eq. flU), all 
the solutions for the mass matrix, and especially the -61,2,3,4 solutions lead naturally 
to such a rather important contribution to dark matter. Inversely the dark matter 
constraint prefers a lower value of the m ee or a higher cos 29 for the solutions -81,2,3,4- 
Note also that for the SMA solution with the central values Sm 2 ol = 5.2 ■ 10~ 5 eV 2 
and tan 2 9 = 1.8 ■ 10~ 3 we obtain rh = m ee to a good approximation which means for 
m ee = 0.39 eV that Q v h 2 = 1.5%. 



Textures and perturbations in the degenerate case: In contrast with solutions ^i, 2,3,4 
the dependence in 9 is already non-vanishing even without any perturbation propor- 
tional to 8m 2 ol and Sm 2 tm . So even in this limit there is in fact a "continuum" of zero 



textures depending on the value of 9, see Eqs. ( pT|) and (0). This allows for a contin- 
uum of basis where the zero texture is diagonal, i.e. of basis which can be considered 
as natural in the symmetry limit. Eventually for some particular values of the angle 
9 the corresponding diagonal mass matrix will display interesting properties. These 
considerations are beyond the scope of this article. 

The structure of the perturbations appear to be also more complicated than for 
^.1,2,3,4- In particular, now there is a dependence of the perturbation from 5m 2 tm in 9, 
which would lead to the interpretation that this perturbation cannot be independent 
of the other perturbations, unlike in the solutions Ai t 2.3,4 where they are independent. 
In the continuum of possible zero textures for the B\ to B4 solutions, a natural choice 
is to take the textures with 9 = or 9 = 7r/4, and make a perturbation around these 
cases (for the SMA and the LMA respectively). In this case the deviation of 9 from 
these values can be considered as a perturbation. For 9 = we obtain: 
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(27) 
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\ 

1 -1 
-1 1 

sin 20 v/2 sin 20 \ 
-2 cos 20 -2 cos 20 
-2 cos 20 -2 cos 20 / 



Here we recover the textures CI and C2 of Altarelli-Feruglio |TTJ(for the single max- 
imal mixing case) and the absolute neutrino mass in the limit Sm 2 ol = 5m 2 atm = 
is m ee / cos 20 ~ m ee . The SMA solution could be obtained from this texture but 
in this already said above, the structure of the perturbation appears to be 

quite complicated, more complicated than from the texture in solutions ^2,3,4. 

We then consider the solution around = ir/4. In this case there are interesting 
zero textures. As explained above, the absolute neutrino mass, in the limit 8m 2 sol = 
0, is now given by m = m ee / cos 20. In this case m ee can be considered as 



5m 



at in 



a perturbation with respect to m: 



( 



M 



B U B 2 



m 



1 

V2 
1 
V2 



1 
V2 
_ 1 
2 

1 
2 



1 

s/2 

1 

2 
_1 
2 
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+ m 
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cos 29 


sin 20-1 


sin 20-1 


\ 




sin 20-1 


cos 20 


cos 20 






V2 


2 


2 




\ 


sin 20-1 


cos 20 


cos 20 


/ 


V2 


2 


2 






/o 





\ 



± rh cos 2 26* 



£-atm 



m cos 26* 



/ 



-1 

V o +i 





/ 



+1 
-1 

-2 cos 29 



\/2sin2# 
y \/2sin20 -2cos2# 



>/2siii20 \ 
-2 cos 29 
-2cos2# / 



(29) 



A* 



B 3 ,B A 



m 



+ m 



/ 75 

1 1 

2 

1 _ 1 

/ cos 29 

sin 20-1 



1 

2 / 

sin 20-1 sin 20-1 



V 



V2 
sin 20-1 



V2 
cos 20 
2 

cos 20 



V2 
cos 20 
2 

cos 20 



V2 



\ 
J 



± 



~ ^atrn 

m — ■ — 



cos 2 29 




1 



+ rh cos 29 



2 



-1 
1 

v 7 ^ sin 20 
-2 cos 29 
-2 cos 29 



\ 
J 



y/2sm 29 
-2 cos 29 
-2 cos 29 



\ 
J 



(30) 



(° 



\o -i 

/ 

\Z2sin20 
\ \/2sin20 

Here we recover the CI and C2 textures of Altarelli-Feruglio (for the bimaximal 
mixing case) and see how the LMA solution can be obtained adding a perturbation 
from these textures which are the only ones where the neutrino masses in the limit 
3 m atm = ^ m2 S oi = doesn't coincide with m ee but can be much larger than m ee . 
Therefore these textures are the most relevant ones for dark matter. In the limit 
of exact bimaximal mixing and m ee fixed this neutrino mass is infinite. However 
as Eqs. (|29|) and (|30|) show, from these two textures we can get the LMA solution 
(with 9 < 7r/4) from the exact bimaximal case considering TTl ee clS db perturbation 
with respect to rh — m ee / cos 29. Because of this involved nature of the solution, 
we could get a LMA solution in this way even if the 11 element vanishes in the zero 
texture, which wouldn't be so easy to see starting directly from the texture zero as 
in the analysis of Ref . [O] . 
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The inverse hierarchy case: We could have an inverse hierarchy solution, only if m ee 
turns out to be not much larger than the lower edge of the range of Eq. ([!]), i.e. 
m ee = 0.05 eV. This is possible only with solutions -82,4, because as with solutions 
A 2t 4,, Sm 2 tm nas to come with an opposite sign with respect to the term in m 2 e in 
m.3, Eqs. (|l7|) -(|T8l). For these solutions -82,4 an inverse hierarchy requires m 2 ~ 



ee 



cos 2 295m 2 atm . For the central values 5m 2 atm = 2.9 ■ 10~ 3 eV 2 § and tan 2 9 = 0.37 §, 
this requires for the LMA solution m ee ~ 0.025 eV. Varying tan 2 6 within the range 
0.2-0.8 || and Sm 2 tm within the range (1.5 — 4.8) • 10~ 3 eV 2 this requires: 



5 ■ 10~ 3 eV < m ee < 5 • 10~ 2 eV (for solutions B 2A and LMA) 



(31) 



For the SMA solution since cos 26 is very close from unity the corresponding allowed 
range coincides with the one of Eq. (|23|). Note finally that if we develop the mass 
matrix around 9 = and 6 = 7r/4, considering m 3 and -1/ 5m 2 ol as perturbations with 
respect to m and a/ Sm 2 tm we obtain the following zero textures: 




\ 

-1/2 
-1/2 / 



and 



/ l/y/2 l/y/2 \ 

1/V2 
V Vy/2 / 



(32) 



respectively (that is to say the textures Bl of Ref. [|TTjl ) • 



7) Ambiguities of the perturbations: In Eqs. (p!9|)-(p2|) we have given the structure 
of the perturbations which has been determined directly by Eqs. (|7|)-(|T8l). They are the 
minimal perturbations which have to be added to give the atmospheric mass difference 
together with the solar angle and solar mass difference. It must be noted, however, that 
these perturbations are not unique. For example in Eq. ( |T9"D we could perfectly add a 
perturbation which is simply proportional to e a tm and the identity matrix. To add this 
perturbation corresponds to a redefinition of the mass scale of the zero texture by a small 
amount. More generally there exist perturbations proportional to e so i we could add which 
are not proportional to the zero texture and which could leave the condition m\ — m\ = 
5m 2 oi unchanged while changing the scale of the zero texture or the atmospheric mass 
difference by a negligible term proportional to Sm 2 ol . To add such a perturbation leads to 
a mass matrix still fulfilling all the experimental constraints, but with another structure 
of perturbation. Note that the atmospheric mass difference is anyway far from being 
determined at the level of accuracy of the solar mass difference. This leads to some freedom 
in the perturbations which exists also for perturbations proportional to 5m 2 atm . We leave 
this discussion to a subsequent publication. Here we gave just the "minimal" perturbations 
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and we do not consider these extra perturbations which are not necessary to satisfy the 
experimental constraints. 

Although the atmospheric neutrino mixing is allowed to be in the range of sin 2 29 a t m = 
0.87 — 1, we considered only the maximal mixing and not the entire range. This and the 
possible small deviations from the constraint U e 3 = may allow for few more perturbations, 
which we disregarded here. 

8) Summary: In summary, instead of starting the analysis from the possible textures, 
we started with the experimental inputs assuming in addition that m ee is larger than the 
neutrino mass differences. We then determined the neutrino mass matrices and in the 
process got the possible textures. We found that, for maximal atmospheric mixing and 
U e3 = and for a given set of values of m ee , Sm 2 ol , 5m 2 atm and 9, there are 8 solutions 
for the mass matrix. Among them the solutions A\ to A± have the interesting property to 
allow a decomposition which is just a zero texture depending only on a single mass scale, 
plus a perturbation depending only on 5m 2 atm , plus a perturbation depending only on Sm 2 ol 
and 9. There is also a set of solution Bi to whose structure of perturbations is more 
complicated but which allow neutrino masses larger than m ee . All solutions and especially 
-Bi, 2,3,4 lead naturally to neutrino masses relevant for dark matter. For solutions -81,2,3,4, a 
lower m ee and a higher cos 29 is preferred by the dark matter fits. We observed also that 
for solutions -81,2,3,4 there is a continuum of possible zero textures. For each solution we 
displayed the allowed ranges for the mass eigenvalues together with the corresponding zero 
texture and we showed how, in the degenerate neutrino mass scenario, the perturbations 
come up to give the full mass matrix. 
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Last minute note : in the process of submitting this manuscript, some of the results derived 



here turn out to have been presented independently in Refs. |22|, |23 [. 
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